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We study the existence of the (thermodynamic) limit of the scaled 
cumulant-generating function L n (z) = |W n | _1 log Eexp{z|Hn W n \} of 
the empirical volume fraction |S fl W n |/|Wn|, where | • | denotes the 
d-dimensional Lebesgue measure. Here E = (J. >1 (Si + X t ) denotes a 
d-dimensional Poisson grain model (also known as a Boolean model) 
defined by a stationary Poisson process n.x = ^\ >;l with intensity 
A > 0 and a sequence of independent copies Si,S 2 ,... of a random 
compact set So. For an increasing family of compact convex sets 
{W„, n > 1} which expand unboundedly in all directions, we prove 
the existence and analyticity of the limit lim^^oo L n ( z ) on some disk 
in the complex plane whenever Eexp{a|So|} < oo for some a > 0. 
Moreover, closely connected with this result, we obtain exponential 
inequalities and the exact asymptotics for the large deviation proba¬ 
bilities of the empirical volume fraction in the sense of Cramer and 
Chernoff. 


1. Introduction and main results. The Poisson grain model (PGM; also 
known as the Boolean model ) is the best studied and most used random set 
model to describe systems of randomly distributed and irregularly shaped 
clumps in a Euclidean space R d , d> 1 [see Matheron (1975), Hall (1988) or 
Stoyan, Kendall and Mecke (1995)]. It is the basic model in stereology and 
stochastic geometry. Statistical analysis of a stationary PGM is mostly based 
on a single realization of the union set of clumps in some region W which 
is assumed to expand unboundedly in all directions [see, e.g., Molchanov 
(1997)]. To be definite in describing our problem, we first give a rigorous 
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definition of a stationary PGM as the union set 

(1.1) S:=(J(S i + X i ) 

i> 1 

of independent copies Si, 52,... (grains) of a random compact set So (typ¬ 
ical grain) that has distribution Q, where the grains are independently 
shifted by the atoms X i, X-j ,... (germ points) of a stationary Poisson pro¬ 
cess n A = X)i>i ^Xi with intensity A (=mean number of germ points in the 
unit cube [0, l) rf ). Throughout this paper, all random elements are defined 
on a common probability space [17,21, P] and E denotes the expectation with 
respect to P. In particular, So is a measurable mapping from [17,21, P] into 
the space of nonvoid compact subsets /C of R d equipped with the Hausdorff 
metric and Q coincides with the image measure P o Sg 1 that acts on the 
corresponding Borel a-field 23(/C) [see Matheron (1975)]. Note that 5 is a 
closed set (P-a.s.) if E|So + B r (o)\ < oo for r > 0, where B r (x ) denotes the 
closed ball with radius r > 0 centered at x S R d and | • | denotes the Lebesgue 
measure in R d [see Heinrich (1992)]. 

The main aim of this paper is to prove the existence and analyticity of 
the limit (as n —> oo) of 

L n (z) := wttt log Eexp{z|S n H7 n |} on D A := (zeC 1 : |z| < 1/A} 

| V* n | 

(1.2) 

for some 0 < A < oo provided that an exponential moment of the volume | So | 
exists, that is, 

(1.3) M(a) := Eexp{a|So|} < oo for some a > 0, 

and (W n ) is a convex averaging sequence of sets in R d , that is, each W n is 
a (deterministic) compact convex set, ( W n ) is nondecreasing and its union 
is R d [see Daley and Vere-Jones (1988)]. Because of the conspicuous anal¬ 
ogy to similar problems in statistical physics [see Ruelle (1969)], we call 
L(z ) = linin^oo L n (z) the thermodynamic limit of (the thermodynamic func¬ 
tion) L n (z). The second aim, which is closely connected with the first, is to 
derive inequalities and asymptotic relationships (in the sense of Cramer and 
Chernoff) for probabilities of large deviations of the empirical volume frac¬ 
tion p n := |5 n W n |/1 W n | from its mean p := E|5 n [0, l) d | = P(o € S). 

In the special case of a bounded typical grain, that is, So C Br{o) for some 
0 < R < oo, both problems were solved satisfactorily by Gotze, Heinrich 
and Hipp (1995) using the device of m-dependent random fields with block 
representation. The proving technique in the present paper is completely 
different from that in Gotze, Heinrich and Hipp (1995) and does not require 
any strong mixing properties of the PGM (1.1) as one would expect. In 
general, (1.3) does not imply specific mixing rates as needed, for example, 
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in Mase (1982) or Heinrich and Molchanov (1999). However, in case of a 
spherical typical grain, (1.3) induces an exponentially decaying /3-mixing 
coefficient [see Heinrich and Molchanov (1999)]. Note that (1.3) does not 
even imply the closedness of S in general; see the Appendix. For a positive 
random variable X with infinite mean, the typical grain Eg = [0, X] x [0,1/A] 
exhibits such an example for d = 2. 

For this reason we choose the probability space [12,21, P] (for its exis¬ 
tence, see the Appendix) in such a way that the mapping R d x 12 3 (x,uj) >—> 
ls(w)0 :c ) i s 23(i? d )<g> 21-measurable. This property allows us to apply Fubini’s 
theorem to the 0-1-valued random field £(x) = l~(x), x E R d , and implies 
that the function 


(1.4) p£\x\, ...,Xk):= E£(xi) • • -£(x fc ) = P(xj E S,.. .,x k G E) 

is 23(i? rffc )-measurable for any k > 1 and |En W\ = f w £(x)dx is a random 
variable over [12,21, P] for any bounded W E 23(i? rf ). The functions (1.4) 
are expressible (and vice versa) by the corresponding probabilities for the 
complement set E c : 


(1.5) 


Pec (xi, .. .,x k ) := E(1 - £(xi)) • • • (1 - £(x k )) 
= P(sn {xi, ...,x k } = 0 ). 


Since (E,; + Xi) n {xi,..., x k } = 0 iff Xi ^ (—S*) + {xi,..., x k }, the shape of 
the probability generating functional (A.2) (of a stationary independently 
marked Poisson process n^g) for v(x,K) = 1 — 1 (~K)+{x 1 ,...,x k ( x ) yields 

Pe°(x i,-..,®fc) = E ITC 1 ~ 1 (-3 i )+{ii,...,* fc }(^i)) 

i> 1 

( 1 . 6 ) 



(k) 

Note that (xi,..., x k ) = 1 — 7s({xi,... ,Xfe}) for an arbitrary random 

closed set E with capacity functional T= [see Matheron (1975)]. The study 
of the sequence (1.2) is closely related with the behavior of the higher-order 
mixed cumulants 

(1.7) d£\xi,...,x k ) :=r(£(xi),...,£(x fe )) for k > 1 

of the random field {£(x), x E R d }, where the mixed cumulant ( semi-invariant ) 
of any random variables Y\,.. . ,Y k (having a finite fcth moment) is defined 
by 


(1.8) T(Y 1 ,...,Y k ):=i 


;-k 


d k 


ds i■■ • ds k 


log E exp < 


12^8 

3 = 1 


■jYj 


ai=.«=Sfc=0 
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and Tfc(y) := T(Y,..., Y) [obtained by putting Y = Y\ = ■ ■ ■ = Y& in (1.8)] 
denotes the kth cumulant of Y. Directly from (1.8) it is seen that, for k> 2, 

(1.9) d£)(x 1 ,...,x k ) := T(1 - £(xi),...,l -£(x k )) = (-l) fc 4^(xi,...,x fc ). 
We are now in a position to formulate our main result. 


Theorem 1. Let E be the PGM (1.1) with compact typical grain Ho 
satisfying (1.3) and let {W n , n > 1} be a convex averaging sequence in R d . 
Then, for any k> 2, 

(1.10) [ \c ( ~\o,x 2 ,...,Xk)\d(x 2 ,...,x k )<(k-l)lH(a)A(a) k - 2 , 

where H(a ) := 8AM(a)(l + exp{AE|Ho|})/a 2 and A (a) := 8(a + A M(a)) x 
(1 + exp{AE|Ho|})/a 2 . Furthermore, the limit L(z) = \im. n ^ 00 L n (z) exists 
and is analytic on the open disk D^ a ) ■ 


The next result states Cramer’s large deviations relationships for the ran¬ 
dom sequence |S n W n \ and an optimal Berry-Esseen bound of the distance 
between F n (x) := P( \J\W n \(p n — p) < xa n ) and the standard normal distri¬ 
bution function <J>(x) = ff oc exp(—t 2 /2)dt/y/2ir, where 

2 _ Var(|SnW n |) 

' |W n | 

= / ^ ——(exp(—A|H 0 U (E 0 — ®)|) — exp(—2A|H 0 |)) dx. 

JR d \W n \ 

The following Theorem 2 is derived from (1.10) combined with a well- 
known lemma on large deviations for a single random variable discussed by 
Statulevicius (1966) [see also Saulis and Statulevicius (1991), Lemma 2.3]. 


Theorem 2. Let the assumptions of Theorem 1 be satisfied and, in ad¬ 
dition, let E|So| >0. Then u 2 converges to a nonzero limit f Rd c£\o,x)dx 
and, for 0 <x< cr nA /| W ra |/2A(a)(l + 4iL n ) with H n = F[(a)/2a 2 l , the asymp¬ 
totic relationships 


n 1_ - F ’n(*) 

(1.11) -- ^ = exp-i 


3 °° 

2T (n) 


1 — $(x) 


[ n 

f k 


X 


1 + 0 


and 


(L12) f4r =exp l 


3 co 

-X ^ (n) 


(n 


—X 


1 + 0 


1 + X 

V^Wni 


1 + X 


VWn 
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hold as n —* oo, where the coefficients 


(1.13) 



1 

( k + 2) (k + 3) 


x 



rfc <+ 2(|sn w n |) 

a n\W n \(ki + 1 )! 


ki> 1, i=l,...,l 

satisfy the estimate \p^\ < 4H n A(a)(2A(a)(l + 4H n )) k /(k + 2)(k + 3) for 
k > 0. Furthermore, there exists some constant c > 0 [ depending on a, X, 
M (a) and <7 2 ] such that 


(1.14) 


sup | F n (x) — <3?(x)| < 
x£ R 1 


c 

viw 


Our next Theorem 3 provides large deviations inequalities for the unbiased 
estimators 

„ . \Enw\ a / \ |sn(H-s)nw r | 

(1.15) Pw :=——— and C w (x) := - — -, x 6 R d , 

of the volume fraction p = P(o G S) and the covariance C(x) = P(o G S, x E 
H), respectively, in the case when the PGM (1.1) is observed on a sampling 
window W £ *B(i? d ). Note that, in contrast to the volume fraction p, the 
covariance C(-) reveals information on the inner structure of the random 
set E [see Matheron (1975) and Stoyan, Kendall and Mecke (1995)]. The 
deviation of the estimators (1.15) from their means p and C(x) is estimated 
under finite-order as well as exponential moment assumptions put on the 
volume of the typical grain Eq. 


Theorem 3. Let E be the PGM (1.1) with compact typical grain Eo 
that satisfy E|Eo| s < oo for some real s > 2. Furthermore, let W C R d be 
a bounded Borel set with inner points. Then there exist positive constants 
ci^(A) and ci^(A) (depending on A and the moments E|Eo| fc , k = 1 ,..., [s], s) 
such that, for any e > 0, 

(1.16) P(| Pw-P\>e)< c( 1 )(A)e- s |lT|- s / 2 
and 

(1.17) P{\C w {x)-C(x)\>£)<cf\X)£- s \W\- s / 2 forallxeR d . 

If Eo satisfies condition (1.3), then the Bernstein-type inequality 

P (pw~P>z) 


(1.18) 


< 


exp 


exp 


1 -P 
'2H(a) 

P 

2A(a) 


■e 2 |lP| 

e\W\ 


if 0 < £ < 
if£> 


H(a)p 


A(a)(l -pY 
H(a)p 


A(a)(l -pY 
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holds for any 0 < p < 1 and H(a),A(a) from Theorem 1. Exactly the same 
hounds hold for the probability P(pw — p < — e). 

In Theorem 4 below, we derive a Chernoff rate function [see Dembo and 
Zeitouni (1998) and references therein] for the sequence of empirical volume 
fractions p n in terms of the thermodynamic limit L(z), which provides an 
extension and refinement of the relationship (1.12) for the x values x(e) = 
e VlWnl/oVi with e £ (0,£*), where e* is determined by the slope of the 
function L(h ) at h = l/A(a). 


Theorem 4. Under the assumptions of Theorem 2 the large deviations 
relationship 


(1.19) 


logP ip n -p>e) 
lirn - —— - 

1^00 \W n \ 


inf 

0</i<l/A(a) 


g(h) = g(h 0 (s)) 


holds in the interval 0 < e < e* := lim^i- 1 /^( a ) L'(h) —p, where g{h) := L(h) — 
h{e + p) and ho(s) is the unique solution of the equation g'{h) = 0, that 
is, L'{h) = e + p. A corresponding relationship is valid for the probability 
P(p n — V < — e), where the function g[h ) is defined for h £ (—1/A(a),0] and 
with —e instead of e. 


This result touches the question of whether p n satisfies the large deviation 
principle, the answer to which seems to be unknown so far. Without giving 
details we mention only that the limit linin^oo L n (h) exists, on the negative 
real axis, which can be shown by the methods of Ruelle [(1969), Chapter 3.4]. 
For related problems concerning large deviation principles for stationary 
independently marked Poisson processes, refer to Georgii and Zessin (1993). 
Similar results for Young measures related to Poisson grain models have 
been proved by Piau (1999). 

The rest of this paper is organized as follows: In Section 2, we investi¬ 
gate (1.2) for a quite general random set model (Lemmas 1 and 2) and put 
together the required tools from point process theory presented in a rather 
general setting (Lemma 3). In Section 3, we are concerned with the proof of 
Theorem 1, which is divided into several steps (Lemmas 4-7), whereas the 
proofs of the Theorems 2, 3 and 4 are deferred to Section 4. The Appendix 
contains, among other things, the construction of a measurable random field 
f(x) = l=(x), x £ R d , and a criterion for (non-)closedness of the PGM E 
given by (1.1). 

2. Preliminary results and relationships to point processes. We first in¬ 
vestigate the behavior of the cumulants Pfc(|E n W n |) and give a condi¬ 
tion which guarantees the existence and analyticity of the limit of (1.2) for 
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the support set Z = supp(£) of an arbitrary (*B(R d ) <g> 2l)-measurable, 0-1- 
valued, stationary random field {£(x),x G R d }. We use the same notation 
as in Section 1. Lemma 2 states that this condition can be expressed by the 
total variation of the reduced cumulant measures of the Cox process 

(2.1) n^^i-i^))^, 

i> 1 

which is directed by the random measure z f R d l(.\(x)(l — f(x)) dx = z\Z c D (-)|, 
where H z = X)z>i is a stationary Poisson process with intensity z> 0 that 
is independent of H. In the second part of this section we introduce a family 
of correlation measures for arbitrary stationary point processes and derive 
(Lemma 3) a recurrence relationship for the corresponding Lebesgue density 
functions provided they exist. Lemma 3 is the key to prove Theorem 1 and 
it seems to be of interest on its own. 


Lemma 1. Let {£(&), x G R d } be a measurable, 0-1 -valued, stationary 
random field on R d with support set H := {x G R d : £(x) = 1}. Then, for any 
bounded W G and k> 2, we have 


\r k (\znw\)\<\w\G k (z) 

with G k (Z) := / \c ( ~\o,X2,...,x k )\d(x 2 ,...,x k ). 

Furthermore, let (W n ) n >i be a convex averaging sequence. If G k (Z) < oo for 
some k> 2, then 


( 2 . 2 ) 


lim 

n —>og 


r fe (|s n w n 
\w n \ 


/ ci fc) (o, x 2 , •.., X k ) d{x 2, ...,x k ), 

J( R d)k-l - 


and if G k {Z) < k\HA k 2 for some H, A > 0 and any k>2, then the thermo¬ 
dynamic limit L(z) = linin^oo L n (z) exists and is an analytic function on the 
open disk Da, where L n {z) is given by ( 1.2 ) with H := {x G R d :f(x) = 1} 
[instead of (1.1)]. For z G Da, the function L(z ) admits the power series 
expansion 


L(z) 




— p 2 ) dx 



(o,x 2 ,...,x k )d(x2,...,x k ), 


where p := p^(o) (volume fraction of Z) and C(x) := p^\o, x) (covariance 
ofZ). 
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Proof. Using Fubini’s theorem and the definition (1.4) we may write 


Eni-n^= E n/ t{xi)d Xi 


(2.3) 


2—1 


i=l" Bi 


= / p^\x 1 ,...,x k )d(xi,...,x k ). 

J B\ x--xB k 

A direct calculation of the logarithmic derivatives in (1.8) leads to 

(2.4) r(y 1 ,...,y fc ) = E(-i) J ’- 1 0'-i)! E iW II 




KiU---LlKj=K i=l \ki£Ki 


[see, e.g., Saulis and Statulevicius (1991)], where the inner sum is taken 
over all decompositions of K = {1,..., k} into j disjoint nonempty subsets 
Ki ,..., Kj. From (2.4) and (2.3) and by repeated application of Fubini’s 
theorem, we see that the integral Js 1 x---xB k r(£(xi), • • • ,C( x k)) d(x i,..., x k ) 
coincides with r(|En B \\,..., |Sn B k \). This means, setting B\ = ■ ■ ■ = B k = 
W and using (1.7), that 

(2.5) r fc (|Sn W\) = [ c^\x 1 ,x 2 ,...,x k )d(xi,X 2 ,...,x k ). 

Jw k “ 

The stationarity of the random field {£(x),x € R d } implies the invariance of 
the mixed cumulants (1.7) under diagonal shifts, that is, 

(2.6) c^\xi,x 2 ,..., x k) = c^\o,x 2 -xu...,x k -xi), 
whence, by substituting yj = Xj — x±, j = 2, ..., k, it follows that 

r fc (|Sn W\) = [ [ c^\o,y 2 ,...,y k )^w(xi) 

J(Rd ) k ~ 1 JR d 

k 

X n 1 i' / (2/l - x 1 )dx 1 d(y 2 ,...,yk) 

1=2 



X \w n (W — a; 2 ) n • • • n (W — x k ) \ d(x 2 ,.. •, x k ), 


proving the first part of Lemma 1. The limit (2.2) is an immediate conse¬ 
quence of Lebesgue’s dominated convergence theorem and the fact that, in 
view of the geometric properties of the W n ’ s [see Fritz (1970)], 

r \w n n (w n - an) n • ■ • n (w n - x fc _i)| 

hm --= 1 

n—>0 ° \W n | 


for any fixed x k -\ E R d . 
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The power series expansion of (1.2) is 
L n (z) =pz + Y, 


r k (\~nw n \)z k 


k> 2 


I W n 


k\ 


and, hence, by our assumptions, 
I L n {z) pz\ <J2 G k(-) 

k> 2 


k\ 


<|z| 2 ffV(|z|A)*~ 2 = , |Z ( g A for|z|<4 

- 11 1 - |z|A A 


Thus, for any n > 1, L n (z) is analytic on the open disk Da and, by (2.2), 
L n (z ) converges to L(z ) uniformly in any closed subset of Da, proving the 
analyticity of L(z) on Da- □ 


To obtain estimates of the form G n (E) < n\H A n 2 in the case of the 
PGM (1.1), we first show that (— z) n c~ l ' > coincide with the nth-order cu- 
mulant density of the Cox process (2.1). In the second step we introduce 

a family of correlation measures 7 ^ n,n ' ) and their Lebesgue densities, and 
study them for 'll = Ilk 1 . In Section 3 we perform a somewhat involved and 
rather lengthy inductive estimation technique to derive bounds of the total 
variation of these correlation measures in terms of moments of |Ho| when 
H is given by (1.1). The basic idea of this method goes back to Ruelle (1964) 
[see also Ruelle (1969), Chapter 4.4], who developed it (without using the 
terminology of point processes) to prove the existence of thermodynamic 
limits for grand canonical Gibbs ensembles with pair interactions. An ex¬ 
tension to ensembles with higher-order interactions was tried by Greenberg 
(1971), but it fails in our situation. 

To begin with, we briefly recall the definition of the nth-order factorial 
moment (and cumulant) measure (and 7 ^) of a point process 'k = 
Y^i>i^Zi that satisfies E'k n (A) < 00 for K e K, by means of its probability 
generating functional = E(]~[j>i w(Z{)), where w:R d 1 —» [0,1] is Borel 

measurable such that l — w has bounded support [see, e.g., Daley and Vere- 
Jones (1988)]. Setting w^;;;f n n {x) = 1 + EjU (vj ~ for 1 - ± < 

Vj < 1, i = 1,..., n, and bounded B 1 ,..., B n 6 ’ ! B(R d ), we define 


a 


(n) 



:= lim 

Vl,...,Vn 


gn 

n dvi--- dv n 




,V n i 



lim dn 

«i,...,unTi dvi ■ ■ ■ dv n 


log Gy 


Bn 


and 
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If af) 1 (resp. 7 ^) is absolutely continuous with respect to the Lebesgue 
measure on R dn , then we denote the corresponding (factorial) moment (resp. 
cumulant density) by p[^ (resp. c^). In the sequel we often write p^\x n ) 

instead p (xq,..., x n ), where X n stands for the (unordered) point set {x\ 

In case the point process is stationary, there exists a unique (signed) mea¬ 
sure 7 v j"| ed on Q3(i? rf ^ n_1 ))—called nth-order reduced cumulant measure — 
such that 

(-Bj-x)) dx 

for any bounded ,B n G *3(R d ). 



Finally, a stationary point process 'h is said to be Brillinger-mixing [see, 
e.g., Ivanoff (1982) or Heinrich and Schmidt (1985)], if E’F ri ([0, l) d ) < 00 and 

the total variation var( 7 ^. ed ) on R d ( n ~ 1 ') is finite for all n > 2 . 


Lemma 2. Let E be the support set of a measurable, 0-1-valued, sta¬ 
tionary random field {£(x),x; G R d }. Then the nth-order reduced cumulant 

measure 7^1) d °f the Cox process (2.1) exists for any n > 2 and its total 

variation (if it exists) takes the form 

var( 7 ^L ) = z n G n ( E) 

V 'n^ ) ,red / V ’ 

with 7^ } (B) = (-z) n [ cf l) (o,X2,...,x n )d(x 2 ,...,x n ) 

n a > red JB 

for B G %$(R d ( n ~~ l ' > ). Consequently, is Brillinger-mixing iff G n (E) < 00 
for all n > 2 . 


Proof. From the shape of the probability generating functional of a 
Cox process directed by an arbitrary random measure [see Daley and Vere- 
Jones (1988), page 262], we deduce that 



which in turn, using the above definitions of moment and cumulant mea¬ 
sures, provides 


a 


( n ) 

nip 



= z"Enis c n^ 

3 = 1 


{n 

nW 



z n r(|E c nH 1 |,...,|E c nH n [). 


and 
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Hence, repeating the steps in the proof of Lemma 1 that lead to (2.3) 
and (2.5) (with 5 C instead of E), we recognize that, for B £ 53(f? dn ), 

a£l ) (B) = z n [ p { ~}(X n )dX n and y ty z) (B) = z n [ (X n ) dX n , 

n 3 JB “ n a JB “ 

where X n = {ay,..., x n } and dX n = d(x i,..., x n ). Thus, the nth-order cu- 

mulant density of ni"'* equals z n ch). The proof is completed by appealing 
to (2.6), (2.7) and the very definition of total variation. □ 


We now introduce a further family of (signed) measures y^"'’ ri on ^8(R rn+n ) 
for n,m > 0 associated with the point process 'L, which is assumed to admit 
moment measures of order m + n. For bounded Hi,..., A m , B i,..., B n G 
*B(R d ) define 

m n 

X AiX X Bj 

i=l j =1 



u i, 

Vi 


Qn+m 

dv\ ■ ■ ■ dv n du\ ■ ■ ■ du r 

—,'UnT 1 


G*[wut, 





• 5 ^n] 
■Vn \ 


For the sake of distinction, let us call the (factorial) correlation 

measure of order (m, n). In case the moment density p( T ™ +n ' ) exists, y^ m,n ^ has 
a Lebesgue density which we call correlation density of order ( m,n). 

Note that Ruelle (1964) apparently first introduced the densities c^ n,n ^ via 
an algebraic method to study cluster properties of the correlation functions 
of classical gases. 

It is evident that 7 is symmetric in the first m as well as in the second 
n components, but not completely symmetric. By logarithmic differentiation 
with respect to v\ in the above definition of 7 ^ we see that 7 ^ coincides 

with y^ 1 ’"' ^ for n > 1 and, thus, c ^ = c^’ ?t ^ for n > 1 provided the 
densities exist. Moreover, for fixed m > 1 and any n > 1, the relationship 
between factorial moment and correlation measures, 


a 


(m+n) 


X AiX X B~ 


( 2 . 8 ) 


„ i=l 


i=i 


= E 

0CJCN 


X AiX X Bi 


a 


i =1 


j&J 


(n—| J|) 


X Bj 
Vie n\j 


holds, where the summation extends over all subsets J of N with |J| ele¬ 
ments. For reasons of consistency, put 7 ^ ,t ’°^(X™ = 1 H i ) = a^ l \x 1 ^ =l A i ) (= 1 
for m = 0) and y^ 0 ’^ (X n - =1 Bj) = 0 for n > 1. 
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To verify (2.8), let us briefly write D v for d n / dv\ ■ ■ ■ dv n and D u for 
d m /dui ■ ■ ■ du n and put 


f(u , v) = and g(v) = 


Then (2.8) is obtained by applying Leibniz’s rule for higher-order derivatives 
of products of functions to the right-hand side of the identity 


D v D u (f(u,v )) = D v (g(v) ■ h(v)) 


with h(v) 


D u (f(u,v )) 

9{v) 


We conclude this section with a recursive representation of the corre¬ 
lation density c{j, n ' n> in terms of the densities j = 1 ,...,n. 

For notational convenience, we omit the superscripts (if confusion is ex¬ 
cluded) and write cq,(X m , Y n ) instead of (xi, ..., x m ,yi ,..., y n ), where 
X m = {x\, ..., x m } and Y n = {yi, ..., y n } are two disjoint sets of distinct 
points in R d . Furthermore, put X , m _ 1 = X m \ {xi} and let |Yj denote the 
cardinality of a finite point set Y c R d - 


Lemma 3. Let T be a point process on R d with strictly positive factorial 
moment densities p[^ for k = 1,... ,m + n>l. Then we have 

(2.9) c *(x m ,Y n )= (-^^^(^.^(yui^.y^y), 

0CYCY n 


where 


( 2 . 10 ) 


K 9 (X m ,Y) 


E (-D |y 

0CVCY 


P*(V U Xm) 

p*(v ui :.!) 


for m, |T| > 1, Y C Y n 


and K x y(X m , 0) = p x y (X m ) /py (X' rn _ ,) form > 1 and Kq,(0,Y n ) = 0 for n > 

1. 


Proof. The relationship (2.8) reads, in terms of densities, as 
(2.11) p*(X m UY n )= c*{X m ,Y)p*(Y n \Y). 

0CYCY n 

Given the moment density functions p^(Y), Y C Y n , with py(0) = 1, there 
exist unique symmetric functions p^(Y), Y C Y n , withp^,(0) = 1 that satisfy 
the equations 

E 4(^\^0=0 for 0/hC Y n . 

0CVCY 


( 2 . 12 ) 
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By means of the functions p\,(Y) we may invert the “convolution equa¬ 
tion” (2.11) by calculating the sum 


E pl(Yn\Y)p^X m UY) 

0CYCY n 

= E PUYn\Y) E MX m ,V) Pq ,(Y\V) 

0CYCY n 0CVCY 

= E ^(X m ,V) E pMYn\Y)p#(Y\V). 

0CVCY„ Y:VCYCY n 


Since, by (2.12), the second sum in the last line vanishes for all proper subsets 
V C Y n , the whole last line is equal to cq,(X m ,Y n ). Using this identity and 
the relationship 


E (-1 )WK*(X m ,V) 

0CVCY 


py(Y U X m ) 
P*(Y U-Cj) 


for 0 c y C h n 


obtained from (2.10) by using the Mobius inversion formula [see Rota (1964)], 
we may proceed with 


Cvji (A m , Y n ) 

= E Pl(Y n \Y)py(YUX^_ 1 ) E (-l) IVI X*(X m ,V) 

0C.YC.Yn 0CVCY 

= E (-1 )WK*(X m ,V) E P**(Yn\Y)p*(Y UX^) 

0CVC Yn Y-.VCYCYn 

= E (-1 )WK*(X m ,V) 

0CVC Y n 

X E PMYn\V\U)py(U uyuiU). 

0C UCY n \V 

Applying again the above derived identity, we see that the second sum in the 
last line equals c^(U U X' m _ 1 ,Y n \ V), proving the asserted relationship (2.9). 
□ 


3. Absolute integrability of the correlation densities of the Cox pro¬ 
cess ni 1} and proof of Theorem 1. Throughout this section we consider the 
factorial moment and correlation densities pq, and merely with respect to 
Cox process Ili 1 ^ defined by (2.1) for the PGM (1.1). For notational ease, we 
indicate this by omitting the subscript 'L at py, cm and Kq,. Our aim is to 
obtain bounds of the integrals f Rdn |c({o},l^)| dY n [= G n (E) by Lemma 2] 
under suitable moment conditions on |Hq|. To do this, however, our inductive 
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proving technique requires us to estimate the integrals J R d n \c(X m ,Y n )\ dY n 
uniformly in X m G R dm for any m > 1. 

Let X m , X' m _ x and Y n be the finite point sets introduced at the end of Sec¬ 
tion 2. Furthermore, for any finite subset Y C R d , put EoOO := Uyev(Eo — y) 
[and thus Eq(Y) = — y)]- For any n > 1 and Y C Y n , define 

(3.1) S(X m ,Y):= (-l) |V/| exp {E( Xl] X' m _ v V)}, 

0CVCY 

where 

E(x-, U, V) := AE|(S 0 - x) D E C 0 (U) n S 0 (F)| for x $ U 
and, for any Y C Y n _i := Y n \ {y n }, 

(3.2) T(y n ,X m ,Y):= £ (-l)l y l exp{—E(xi,y n ; X' m _ x , V)}, 

0CVCY 

where E(x,y;U,V) := AE|(S 0 - x) n (E 0 - y) n Eg(f7) n E 0 (y)| for x^y, 
xtfiU and y ^ V. Note that E(x, y; U, 0) = E{x\ U, 0) = 0, implying S(X m , 0) 
T(y n ,X m ,0) = 1. 

From (1.6) it is clear that 

p{v\jx" 1 \) = ex P{ -AE l H ol + AE K S o - ®i) n E 0 (V U X^_i)|} 

= exp{-AE[(E:o - Xl ) nE^(X' m _ 1 )\}exp{E(xi-X' m _ 1 ,V)}, 
so that, by (2.10) and (3.1), 

(3.3) K(X mi Y) = exp{-XE\(E 0 -x 1 )nE c 0 (X' m _ 1 )\}S(X m ,Y). 

Next we establish a recursive representation of S(X m . Y n ) with respect to Y n 
in combination with the nonnegative terms T(y n , X m , Y ) for Y C Y n -\. It 
turns out (see Lemma 5 below) that the integrals f Rdn T(y n ,X m ,Y n _i) dY n 
can be represented as functionals of certain PGM (3.6), which enables us to 
derive upper bounds of them under reasonable moment conditions on the 
volume of the typical grain Ho- By means of these bounds and the following 
Lemma 4 we find corresponding bounds of J Rdn \S(X m , Y n )\dY n which in 
turn, using (2.9) with (3.3), enable us to establish the desired bounds of 
f R d„ \c(X m ,Y n )\ dY n ; see Lemma 7 below. 

Lemma 4. We have 

S(X m ,Y n ) = S(X m ,Y n _ i)(l - exp{E(xi;^„ 1 ,{y n })}) 
-exp{E(xr,X' m _ 1 ,{y n })} 

X 5Z T (yn,X m ,Y) 

0CYCY n _ 1 


xexp{E(x 1 ;X^ l _ 1 ,Y)}S(X m UY,Y n _ 1 \Y). 
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Proof. By the definition of the terms E{x\ U, V) and E(x, y\ U, V), and 
the relationship |^4| + \B\ — \AdB\ = \AuB\ for bounded A,B£ we 

get 

E(xi', X^—i, Y U {j/n}) 

= E(x 1 \X' m _ 1 ,{y n }) + E{x 1 \X' m _ 1 ,Y) - E{x 1 ,y n ]X' m _ l ,Y) 
for any Y C Y n -\. Furthermore, we may rewrite the sum S(X m ,Y n ) as 
E ( — l) |5| ( ex p{F ; (xi;X^_ 1 ,y)} — exp{E{x\\X' m _ 1 ,Y U{y n })}). 

0C7C7„_ 1 

This combined with the foregoing relationship leads to 

S(X m ,Y n ) = S{X m ,Y n _ 1 ) 

- exp{Fl(xi; X' m _ 1: {y n })} 

x E (-l) |F| exp {E{x 1 -,X' m _ 1 ,Y)-E{x 1 ,y n -,X , m _ 1 ,Y)}. 

0CYCY n -! 

A simple application of the Mobius inversion formula [see Rota (1964)] to 
the terms (3.2) yields 

exp{-E(xi,y n ;X^_ 1 ,Y)} = l+ E (- 1 ) lul T(y n ,X m ,U), 

0CUCY 

Inserting this identity on the right-hand side of the previous equality we 
arrive at 

S(X m ,Y n ) 

= S(X m ,Y n _ 1 )(l-exp{E(x 1] X’ rn _ 1 ,{y n })}) 

-exp {E(x 1 ;X' m _ 1 ,{y n })} 

x E E (-l) |y| - |t7| T( 2 / ri ,X m ,C/)exp{F;(x 1 ;^_ 1 ,y)}. 

0 c 7 cy „^ 1 0cucy 

By interchanging the sums and substituting V = Y \ U we obtain that 
E E (-l) lYHUl T(y n ,X m ,U)exp{E(x 1 -,X , rn _ ll Y)} 

0CYCY n -! 0CUCY 

= E E (-l) IYHUI T(y n ,X m ,U) 

0CUCY n _ 1 Y:UCYCY n _ 1 

xexp{E(x 1 ;X^ l _ 1 ,(Y\U)UB} 

= E T(y n ,X m ,U) E (—l) 1 ^ 1 exp{E(x\-, U 17)}. 

0 cucY n . 1 0 <zy<zY n - 1 \u 
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Since \A U B\ = \B\ + |D B c \ for any bounded A,B € %$(R d ), 
E(x 1 -X’ m _ 1 ,VUU) = E(x 1 ;X’ m _ 1 ,U) + E(x 1 ;UU X ^_,, V ), 
whence, by definition (3.1), it follows that 

£ (—l) 1 ^ 1 exp{E(xi;X^ n - 1 ,V U U)} 

0CVCY n ^\U 

= exp{E(xi; U)}S(X m U U, Y n _ 1 \ U). 

Finally, assembling all the above identities we obtain the assertion of Lemma 4. 

□ 


Lemma 5. Let E be the PGM (1.1) with compact typical grain Eo satis¬ 
fying E|Eo| n+1 < oo for some fixed n> 2. Then, for any m> 1, 

sup / T(y n , X m ,Y n _i) d\ n 
x m J R dn 

(3.4) < (n — 1)! 

xY"— V E|E 0 | ni+2 E|E 0 | n2+1 E|E 0 [^ +1 

X hi kl n 1+ n 2 +hn k =n-l ^ ^ "* ! 


If condition (1.3) is satisfied, then the estimate 

(3.5) sup/ T(y n ,X m ,Y n _i)dY n <n\(-') 
x m JR dn \aj 

holds for all n > 2 and m> 1. 


2\ n XM(a) ( 1 | A M(a) 


n —2 


Proof. According to the definition (3.2), 

T(y n ,X m ,Y n _ 1 ) = Y, (—1) |Y| Gxp{—E(x 1 ,y n ]X' m _i^Y)}. 

0CYCY n -! 

Now, for any nonempty Y C Y n _ i, we introduce a new PGM E(xi ,y n ; X' m _ 1 ,Y) 
governed by II A = X)j>i and the typical grain (E 0 - an) n (E 0 - y n ) H 

E o( x ' m -i) n EoOO, that is ; 

^(^1: Vni X m _i , Y) 

(3.6) 

:= (J ((Ej - X!) n (Ej - y n ) n n 5 i(Y) + Xi), 

i> 1 

where E i(Y) = (j yeY ( E i ~ v) and ’E c i (X' m _ 1 ) = f|)= 2 ( s f - ®i)- 
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Obviously, for each realization of the PGM (3.6), we have 

E(x 1 ,y n -,X' m _ 1 , Y )= (J E(x 1 ,y n -X' m _ 1 ,{y}) for Y C Y n _\. 

y& 

Applying the well-known formula P(o 6 H) = 1 - exp{—AE|H q|} [which is 
valid for the PGM (1.1)] to the stationary PGM (3.6) we see that 

exp{-E(xi,y n \X’^^Y)} = P(o £ E(xi, y n \Y)) 

= 1-P ( IJ {oGS(*i,i/ n ;X^_i,{y})} 

\ y eY 


Since J2e)CYCY n -1 ( — 1)^ = 0, it follows from the inclusion-exclusion prin¬ 
ciple that 

T(y n ,X m ,Y n _ 1 ) = - (-l) |r| P U{° eS (h,!/»;^-i,M)} 

0CYCY n -i \yeY J 


fn— 1 


n— 1 


P (^n -*-S(xi,yn;V^_ 1 ,{j/ i })(°) J- 
Thus, by Fubini’s theorem, 

j Rdn T (yn,Xm,Y n -l)dY n = t(yj Rd 1 S( I i I | h ;^ 1 ,{y})(o) d y 
and, for each realization of (3.6), 

J Rd -*-=({zi,.M)(°) d V — J Rd ^2 ^-(pi-xi)n(Si-y n )n(Si-y)(^i) dy 


dy n 


i> 1 


— I '“‘® I —xi)r\(E.i—y n ) (-^*)) 


i> 1 


whence, by applying the polynomial formula and using Fubini’s theorem 
again, we get that 

n—1 


'Rd 


*-3(si,J/ n J-X'm-iiM) dV 


< J_ y- ( n 1)- 
“t 1 fc' ril+ ...^ =n _ 1 n 1 !---n A: ! 
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n-1 \k 


X X] n(l(^-,i)n( 3 ,-,n)(^)l^ 

\h,— ,*fc>i j=i 
k 


— (n_l)!^— ^ j-j- E|(r,Q x i) n (-=.q y n )||r,o| n - 7 


k I 

fc=l ' niH-|-nfc=n—lj=l 


njl 


Here the sum Y* stretches over fc-tuples of pairwise distinct indices and 
the last equality is obtained by applying the Campbell-type formula (A.3) for 
fj(x,K) = l(K- Xl )n(K- yri ){x)\K\ n i . Together with the obvious relationship 

(3.7) [ \(Z 0 -x 1 )n(Z 0 -y n )\dy n = \Z 0 \ 2 , 

JR d 

we finally arrive at the desired estimate (3.4). 

The existence of the exponential moment M(a) of |Ho| implies E|So| fc < 
k\M(a)a~ k for all k > 1. Inserting this moment bound in the right-hand side 
of (3.4) and taking into account 


ni-l- \-n k =n —1 


2=1 


V (n 1+ 2)n(n. + l)<C Jb’V, 


we obtain that 


J Rd ^ (yJ Rd (°) dy'j dy. 


< rd y, 1 (XM(a)) k /n- 2\ n 

~ a n fifi a k k\ V k - 1 / 
k =1 


< nil — 


1 AM (a) 


1 + 


AM(a) 


n—2 


This is exactly the desired estimate (3.5). Thus, Lemma 5 is completely 
proved. 


□ 


Lemma 6. Let S be the PGM (1.1) with compact typical grain So that 
satisfies E|So| n+1 < oo for some fixed n> 1. Then, for any m > 1, 

(3.8) sup [ \S(X m ,Y n )\dY n <c n ( A) < oo, 

A' m JRdn 

where the constant c n ( A) depends on A and the first n + 1 moments of |Ho|. 
Moreover, if condition (1.3) is satisfied, then (3.8) holds with 

(3.9) c n (X)=n\A n B(l + B) n ~ 1 
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for all n > 1 and m > 1 with A = |(1 + exp{AE|Ho|}) = and B = 

XM(a) 


Proof. In view of the obvious inequalities 

E ( x A x 'm-i,{yi}) < E|(H 0 - aq) n (H 0 - yi)| < AE|H 0 | 

and e x — 1 < xe x for x > 0 together with (3.7) we see that 

\S(X m ,{y 1 })\dyi= f (exp{p;(aq; X^, {yi})} - 1) dyi 
J R d 

< Aexp{AE|S 0 |} f E|(“ 0 - aq) n (H 0 - yi)| dyi 

JR d 

A 


IR d 


1 ~P 


FI” |2 
-t —n • 


Dehne 


Am,n ■ sup / \S( K X mi Y n )\dY n and B mn . sup J T( y y ni X rni Y n _\) dY r 
Xm J R dn X m J R dn 

for m,n > 1. Prom Lemma 4 and E(x\; Xf l _ 1 , Y) < AE|Ho| for Y C Y n _\ we 
get the inequality 


(3.10) 


Am,n if 


+ 


A 


1 — p 


■El Sr 


n— 1 

E 




n — 1 


ddm,k+l-Am+k,n—k—l 


with A m fi = 1 and A m: i < AE|Ho| 2 /(l — p) for any m > 1. Since, by Lemma 5, 
Bm,k +1 h\Ck for k > 1, where depends on A and the first L + 2 moments 
of |Ho| but not on m, we recognize by induction on n, that -A m>n < n\D n , 
where Dq = 1, D\ = AE|Hq| 2 /(1 — p) and 


n —1 


D n = 


(1 — p)n 


Du-iCq + 




^ ^ CkDn—k —1 


for n > 2 with Co := AE|Ho| 2 . Therefore, A mjn does not depend on m and is 
bounded by terms that involve merely A and E|Ho| fc , k = 1,2,..., n + 1. This 
proves the first part of Lemma 6. 

We also prove (3.9) by induction on n. Prom (1.3) we get E|Ho| 2 < 2M(a)/a 2 , 
implying 

2AM(a) 


Am 1 < < Aexp{AE|“o|}E|H 0 | 2 < 


a 2 (l — p) ’ 
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which is even slightly stronger than (3.9) for n = 1. Assume now the validity 
of (3.9) for n = 1,..., N — 1. Taking into account the estimates B m ^ + 1 < 
k\Ck with Ck = (k + \)B(^) k+l (l + B) k_1 for k > 1 as stated in Lemma 5 
together with Co < 2B/a, we may write 

Dn -N(1-p) 

x ( 'A^Bil + B) n ~ 2 ^ 

1 /o\fc+l \ 

+- V (k + l)B[ - ) (1 + B) k - 1 A N ~ k ~ 1 {l + B) N - k - 1 . 

After a short calculation using that Jfk>i(k + 1)(^ ) A ' = (l-p)(3-p), we 
arrive at 

T>tv< A n B( 1 + B) n ~ 1 for N >2. 

Thus, the second part of Lemma 6 is proved. □ 

Lemma 7. Let E be the PGM (1.1) with compact typical grain So that 
satisfies E|So| n+1 < oo for some fixed n > 1. Then, for any m > 1, 

(3.11) sup / \c(X m , Y n ) | dY n fC c m ^ n (A) <C oo, 

Xm JR dn 

where the constant c m:n (X ) depends on m, X and the first n + 1 moments 
o/|So|. If condition (1.3) is satisfied, then (3.11) holds with 

(3.12) c m , n ( A) = n\2 m+1 AB(4A(l + B)) n ~ l 
for all n > 1 and m> 1 with A and B as in Lemma 6. 

Proof. Replacing K(X m ,Y) in (2.9) with (1.6) leads to 
c(X m ,Y n ) = exp{-AE|(H 0 - a*) n ^(X^)]} 

x y (-i)\ Y \s(x m ,Y)c(Yux l rn _ 1 ,Y n \Y). 

0CYCY n 

Since c(X m , 0) = exp{—A|So(Ai m )|} < 1 by (1.6), c(0, Y n ) = 0 for n > 1 by 
definition and S(X m ,0) = 1 for m > 1, and since both S(X,Y) and c(X, Y) 
are symmetric in Y C Y n for fixed X, we deduce from the latter recurrence 
relationship the inequality 

[ \c(X m ,Y n )\dY n 

JRdn 
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< 


(3.13) 


[ \c(X' m _ 1 ,Y n )\d,Y n + f \S(X m ,Y n )\dY r 

J R dn J R dn 


I R d 
n— 1 


+ Y(l) J^\s(x m ,n)\ d Y k 

x sup [ I c(Y k U X' m _ v Y n \ Y k ) I d(y n \ Y k ). 

Y k JR d ( n ~ k ) 

For any m > 1 we have 

[ \c(X m ,{y})\dy= f exp(-AE|H 0 (X m ) U (H 0 - y)|) 

x (1 -exp(-AE|H 0 (X m ) n (H 0 - y)\))dy 
< A(1 -p)E|S 0 (X m )||“ 0 | < AmE|“ 0 | 2 - 


Using the estimate (3.8) of Lemma 6 and applying (3.13) successively 
to the remaining integrals on the right-hand side of (3.13), we obtain a 
bound of the left-hand side of (3.13) in terms of Cfc(A), k = and 

supy J R d\c(X U Y. {y n })\ dy ni Y C Y n _i, X C X' m _ 1 . This combined with 
the foregoing inequality proves (3.11). 

We now assume (1.3), which gives E|So| 2 < 2M(a)/a 2 , so that together 
with m < 2 m_1 , 


f \c(X m ,{y})\dy<2 m ~ 1 AB, 

JR d 


which implies (3.12) for n = 1 and m > 1. Let now (3.12) hold for all m,n> 1 
that satisfy m + n< M + N. Then, making use of estimate (3.9) of Lemma 6, 
it follows from (3.13) that 



c(X M ,Y N )\dY N 


< N\2 m ~ 2 (4:A) n B(1 + B ) 7V " 1 + N\A n B( 1 + B) n_1 

N -1 

+ N\J2 AkB ( l + B) k ~ 1 2 M+k ~ 2 (4:A(l + B)) N ~ k 


k =1 


= N\2 M ~ l (AA) N B(l + B) n ~ x 


1 | 1 
2 2 2N + m ~ 1 


N—l 


+ E 


i 

2 k + 1 


Thus, the validity of (3.12) for m + n = M + N follows because the sum in 
brackets does not exceed one for M + N > 2. This completes the proof of 
Lemma 7. □ 






22 


L. HEINRICH 


Proof of Theorem 1. As an immediate consequence of (3.12) for 
m = 1 and Lemma 2 applied to the stationary PGM (1.1), we obtain that 

(3.14) G n+1 (~)= [ \c{{o},Y n )\dY n <c ltn {\)<n\AAB(AA{l + B)) n - 1 

J R dn 

for all n > 1. Thus, by the definition of A and B in Lemma 5, we get (1.10) 
with H(a ) = 4 AB and A (a) = 4A(1 + B). Finally, the existence and ana- 
lyticity of the thermodynamic limit L{z ) of the function (1.2) on the disk 
Ha( a ) follows from the second part of Lemma 1. □ 

4. Proofs of Theorems 2 4. 

Proof of Theorem 2. From Lemma 1 and (3.14) we get the estimate 

(4.1) \T k (\EnW\)\<\W\(k-l)\H{a)A{a) k - 2 for k > 2 

and any W £ 03 (R d ). For the standardized random variable := (|Sn W n | — 
p\w n \)/a n vm (with a n > 0 ), (4.1) implies that 

(4.2) |r*(6»)| < (k - l)!^Sp 2 < for k >3 

with H n = H[a)/2u\ and A n = cr„y / |lL 7 n|/A(a). Note that the asymptotic 
variance linin^oo cr^ = J Rd c ^ (o, x) dx is finite and strictly positive iff 0 < 
E| So | 2 < oo. In this case we can find suitable upper and lower bounds of 
ch’(o,x) = exp{—AE|5 0 U (So — a;)|} — exp{—2AE|S 0 |} that lead to the in¬ 
clusion 

exp{—2AE|S 0 |}(1 - exp{-AE|S 0 |}) ^(Z!^ 

tpo| 

< [ c~\o,x)dx < AE|S 0 | 2 exp{—AE|S 0 |}. 

J R d “ 

The estimate (4.2) enables us to apply to ^ n a well-known lemma on large 
deviations of a single random variable proved by Statulevicius (1966) which 
immediately provides the asymptotic relationships (1.11) and (1.12) as well 
as the Berry-Esseen bound (1.14) stated in Theorem 2. To be precise, accord¬ 
ing to the result by Statulevicius (1966), the relationships (1.11) and (1.12) 
are only valid in the narrower interval 0 < x < <5*A n for any 5* < 5q(1 + 5q)/2, 
where 5o £ (0,1) denotes the unique real root of (1 — <5) 3 = 6H n 5. Indeed, 
since H n > 1/2, by (3.14) for n = 1, we have <5o(l + 5q) < 5q/(1 — <5o) 3 = 
1/6 H n < 1/(1 + 4 H n ). Using again (4.1) and the inequality ( fc+ | +1 ) < 2 k+l , 
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we can estimate the coefficients (1.13) as 


idr’i 


< 


fc+i 

£2 


k+l 


(k + 2){k + 3) l=1 


£ n 

.-\-ki=k+l i —1 


fci+- 


A (a) ki H(a) 


(J* 


2 k A {a) k+1 k \(2H(a)\ l 

(k+2)(k+3)^\i- iA ^ ; 


4g w A(g) 

(A: + 2) (A + 3) 


(2A(a)(l+4^)) fc 


for A: > 0. Therefore, the series J2k >o mY i x / a n^/\W n \ ) k converges absolutely 
for |x| < A n /2(1 + 4 H n ), and the O-terms in (1.11) and (1.12) can be easily 
verified by evaluating the remainder terms given by Statulevicius (1966). 
Thus, (1.11) and (1.12) are valid for the whole interval 0 < x < A n /2(1 + 
4:H n ), which completes the proof of Theorem 2. □ 


The proof of the large deviations inequalities stated in the Theorem 3 
relies on Chebychev’s inequality combined with Lemma 7 and (3.14) [resp. 
Lemma 1 and (1.10)]. 


Proof of Theorem 3. For any integer N > 2, the JVth moment of a 
random variable Y can be expressed by its cumulants Tfc(y), k = 1,... ,N 
[by inverting (2.4)] in the manner 


(4.3) 


E Y n 


N 


£ 


N\ 

~k\ 


£ 

niH- \-n k =N 


r ni (Y) 

n\\ 


r» k (Y) 
n k \ 


Consider (4.3) for Y = |Hn W\ — p\W\ with p = E|Sn [0, l) rf |. Since Ti(y) = 
E Y = 0 and, by Lemma 1 combined with (3.14), |T n (y)| < ci jn _i(A)|LF| for 
n = 2,..., N, we are led to 


E Y n 


[A/2] 

<N\ •£, 

k =1 


\wt V 

k\ ^ 

niH- \-n k =N 


c l,ni —l(A) 

n\\ 


c l,n fc -l (A) 

n k \ 


<cY(A)max{|lT|,|lF|[ 7V / 2 ]}, 

where c'^) (A) depends on the first N moments of |Ho|- Hence, for an even 
integer s > 2, (1.16) follows from Chebyshev’s inequality. To prove (1.16) for 
any real s > 2 we next show 

(4.4) E\\EnW\-p\W\\ s <c^(X)\W\ s / 2 for s > 2 
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provided that \W\ > 1. For this, we introduce the “truncated” stationary 
PGM Ew ■= Ui>i(Si V + Xi) generated by IIa = J2i>i^x t and the typical 
grain 

"Jf :=S 0 n B r w(o), 

where the random variable Rq V := sup{r > 0: |Eo n B r {o)\ < |EF| a } takes the 
value oo if |Eo| < \W\ a . Here, we put a = N/2(N + 2 — s) i£N<s<N + 2 
for some even integer N >2. Define Y\\r := \E\y H W\ — |VF|E|E^ fl [0, l) d | 
and let E\y denote the PGM with typical grain Ho \ E$ . Since E\y C E and 
E \ E\y C E\y we have 

\Y\ < \Y W \ + max{|Hiy n W\,\W\E\E w n [0, l) d \}, 


which implies 

(4.5) E|F| S < 2 S ~ 1 E\Y W \ S + 2 s " 1 E|H M / n VF| S + 2 S - 1 (AE|H 0 \ H^||PF|) S . 

By definition of Eq V , 

E|H 0 \ E^\ k = E|H 0 n B c r w ( 0 )| fc l { | Ho |>|^| a} < E|H 0 nVF|-“( s - fc ) 

for 0 < k < s. 


Thus, E|E 0 \E^||W'| < E|H 0 | S |VF| 1_Q: ^ _1 ^ < E|H 0 | s v^. Next, applying (A.3), 

e|h w n it| s < \w\- n e ^ ^((h* \ S V) + * i ) n w\^j 


N 


< |VF| s_Ar lV! ^2 


(\\W\) h 


k =1 


k\ 


e n 


E|H 0 \ Eff \ ni 


niH- \-nk=N i =1 


nil 


< ci(iV)|VlT| s - 7V m^{(AE|3 0 | s ) fc |W'| fc(1 -“ a)+ “ JV } 

<c 2 {N,\)\W\ s / 2 . 


Since, by Lyapunov’s inequality, E|Yve| s < (E|lV| Ar+2 ) s /( Ar+2 ), we need only 
to verify that E\Yw\ N+2 < cn{X)\W\X +2 Y 2 : which in turn follows from (4.3) 
(with Y\,y and N + 2 instead Y and N) whenever |rAr + 2 (|Hyv H W\)\ < 
cs(N, A)|W|(^ +2 )/ 2 . A thorough examination of the proofs of Lemmas 5- 
7 reveals that the constant ci jn (A) in (3.14) takes on the form 


ci,n(A) = AE|H 0 r +1 + 6«( A)E|H 0 | n + b^( A), 


where bn\ A) and bn\ A) are given polynomials in exp{AE|Ho|} and the first 
n — 1 moments of |Eq|. Hence, by E|H( ) f/ | Ar+2 < |VF| Q ( A+2_s ^E|Ho| s , we get 
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the desired upper bound of |r/v+2(|—vp hi W|)|. Putting together the above 
estimates yields (4.4) and hence (1.16) is proved for any real s> 2. 

To establish the second inequality (1.17), we use that H U (H — x) is also 
a stationary PGM with typical grain Ho U (Ho — x) and volume fraction 
p(x) := P(o G H U (H — x)). In view of the obvious decomposition 


C w (x)-C(x) =Pw~P+ 


|(H-x)nlP| 

WT 


-p- 


|(H U (H — x)) n w\ 

\w\ 


-p(x] 


we obtain (1.17) by applying (1.16) to the three stationary PGMs H, H — x 
and HU(H-r). Finally, to prove the exponential inequality (1.18), we again 
employ a Chebyshev-type inequality. In this way we obtain, for e > 0 and 
0 <h< p/ A(a), that 


P(pw ~ P > e) < exp{— h\W\(e + p) + log 


<e W {-h\W\e + ^^\W\^2(hA(a)) k - 2 
l k> 2 


< exp j-h|IT|e + 

Taking h = e(l — p)/H(a) for 0 < e < H(a)p/A(a)(l — p) proves the first 
part of (1.18), whereas the second part is obtained by setting h = p/A(a) in 
the latter inequality. □ 


Proof of Theorem 4. As in the proof of Theorem 3, using (1.2) and 
the notation p n = p\y n , 

P (Pn~P>£) < exp{|W n |(L n (h) - h(e + p))} 

for any h> 0, whence, by virtue of Theorem 1, it follows that 

lim sup 77 - 7 —r log P (p n — p>e) <g{h ) for 0 <h< 

n— >00 \W n \ A (a) 

Thus, the limit on the left-hand side is bounded from above by inf 0 <ft<i/A(a) 9(h)- 
Relationship (1.19) is proved as soon as we show that 

(4.6) liminf ——— log P(p n — p > e) > inf q(h). 

v ' 00 \W n \ K ~ 0<h<l/A(a) y ’ 

For brevity put Cn(e) := |H n W n \ — (e + p)\W n \. Then, for any 5 > 0 and 
h> 0, 

p (Pn-p>e) 


(4.7) > P(p n -p€ (e,e + <5]) 

> exp{|R 7 ri |(L n (/i) -h(e + p + 5))} 


Eexp{hCn(g)}l{c„(£)g(OAlVFn|]} 

Eexp{hCn(g)} 
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Due to the properties of cumulant-generating functions [see, e.g., Dernbo 
and Zeitouni (1998), page 27], the functions L n (-) are convex on the whole 
real axis and L"(h) > 0 for every h G R 1 (provided that > 0). In view 
of Theorem 1, L' n (h) —■> L'{h) and both L' n (h ) (for n > no) and L'(h ) 
are strictly increasing for 0 < h < 1/A (a). Likewise, L'^{h) —* L"(h ) with 

L"{h) > 0 for 9 < h < 1/A(a). Hence, for each e G [0,e*) and sufficiently large 
n, there exists a unique h n = h n (e) G [0,1/A(a)) that satisfies the equation 
L' n (h n ) = e+p. Moreover, we have h n —* ho, where ho = ho{e) is the unique 

solution of L'{h) = e + p. Since h i—> g{h) is a convex function and g'(ho) = 0, 
it follows that g(ho) = inf Q <^ < iM( a ) g(h). Consequently, putting h = h n on 
the right-hand side of (4.7) and taking into account that L n (h n ) —> L(ho), 

n—>oo 

we arrive at 

liminf log P( T ra ~-P - £ ) 
n->oo \W n \ 

(4.8) 

>g{h 0 ) - /i 0 <5 + limmf log(G re (5|W n |) - G n (0)), 

where the distribution function G n (x) = Eexp{h n C n (e)}l^ n ^ <x y/Eexp{h n C n (e)} 
possesses the Fourier-Stieltjes transform G n (t) = Eexp{(*t + /i n )Cn(e)}/Eexp{/i n Cn(£)}- 
Using (1.2) and L' n [h n ) = e + p we can write 


log G n (t) = | W n \(-it(e + p) + L n (it + h n ) - L n {h n )) 

= ~\W n \t 2 ( (1 - i &)Ln(i'&t + h n )d'& if \t\ + h n < 1/A(a), 

Jo 

where the last line is obtained by partial integration of L"(*M + /i n ) with re¬ 
spect to d G [0,1]. An application of Theorem 1 shows that log G n {t/ J\W n \) —» — t 2 x 

n —xx) 

L"(ho)/2 for all t G R 1 , which in turn implies G n (xJ\W n \) —> &(x/WL"(ho )) 

77/ XX 

provided that L"(ho) > 0. In this case, G n (8\W n \) — G n ( 0) —> 1/2, prov- 

TL XX 

ing (4.6) and, thus, the desired relationship (1.19) holds. If L"(/io(eo)) = 0 
for certain eo G (0, e*), then there exists some g > 0 such that L"{ho{e )) > 0 
for e G [so — hi £ o + h\ \ {eo}- Since log P (p n — p > e) is nonincreasing in e, it 
follows that 


g(ho(eo + h)) < lim inf 

n—>oo 

< lim sup 

77 ,—>00 


logP(Pn -V> £p) 
\W n \ 

logP(p n ~p>£o) 

\W n \ 


<9{ho{e o-rj)). 
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Hence, having in mind the continuity of ho(-), the proof of Theorem 4 is fin¬ 
ished. 

□ 


APPENDIX 

For any random closed set E (defined as (21, ctj)- measurable mapping 
S:[fl,2l, P] i—> [tF, af], where T is the (metrizable) space of all closed set 
in R d and cry is its Matheron cr-field), the mapping ( x,lu ) i—► 1 h(^)(®) is 
(23 (R d ) (g> 21)-measurable since (x,F) i—» 1 F (x) is (23 (R d ) <g> (Jj) -measurable 
[see Matheron (1975), Chapter 2]. As we see below, the PGM (1.1) is no 
longer P-a.s. closed in R d if E|Eo| < oo and E|So + B e (o)\ = oo for any e > 0. 

To preserve the (23(i? d ) <g> 21)-measurability of the indicator function l 3 ( w ) (x) 
(needed to apply Fubini’s theorem), we define the not necessarily closed 
PGM (1.1) as a set-valued, measurable mapping of an independently marked 
Poisson process Ha,q := )Ci>i <5 [y t s, ; ] with mark distribution Q(£) = P(Ho E 
£), £e23(/C). 

More precisely, let M£ denote the space of all integer-valued measures 
ip on [R d x /C, 23(if d ) <g> 23 (/C)] that satisfy ip(B x 1C) < oo for each bounded 
B E *3(R d ) and let Tl/c be the cr-field generated by the sets {ip E Mjc : ip(B x 
C) = n} for n > 0, bounded B E 23 (R d ) and £ E 23 (1C). Each ip E M% ad¬ 
mits a representation ip = as a sum °f Dirac measures with 

respect to the at most countable set of atoms [xi(ip),ki(ip)\, i > 1, where 
each atom is counted according to its multiplicity. Note that the mappings 
Mfc 3 ip i—► [xi(ip),ki(ip)] E R d x /C are measurable [see Matthes, Kerstan and 
Mecke (1978)]. Finally, define 

(A.l) Z(ip):=\J(xi(ip)+ki(ip)). 

i> 1 


Proposition 1. The mapping ( x,ip ) e-> £(x,^>) := l 3 (^)(x) is ( f B(R d )^ 

9Jbc) -measurable, that is, {(x, ip) R d x Mjc : x ^ E(?/>)} E 23(i? d ) <g> 

Proof. Let B be the countable set of open balls in that have rational 
radii and midpoints with rational coordinates. For any sequence {K n , n > 1} 
in /C that satisfies K n f R d put S n ( ip ) := Uj:o; i eA' n ( :c i(V’) + ^j(V’))- Obviously, 

W G 1C and E(V>) = U n >i s n(V0- 

It can be readily checked that the set £ -1 ({0}) = {(x, ip) E R d x Mjc :x£ E(ip)} 
coincides with f)n>i Ubgb(B x {V ; G M/c : H n (^)nB = 0}). However, this set 
belongs to 23 (R d ) X DJI K if {ip : E n (ip)nB = 0} = {ip: J2i>i 1 'R. B n(K n xK)(^i(' l P), h(ip)) 
0} E Wl/c for any n> 1 and B E B, where 1Zb := {(x, K) € R d x IC:(x + K) n 
B / 0}. Since the mapping ip i—> f{ x i{' l P)i h(ip)) is 9Jt/c-measurable 
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whenever /: R d x /C i—> R 1 is (53 (R d ) <g) 53(/C))-measurable [see Matthes, Ker- 
stan and Mecke (1978)], we need only to verify that 7 Zb G IB (R d ) (g>53(/C) 
for B £ B. Since there exists a sequence of closed balls B n such that B n } B 
and thus 1Zb u T Rb, it suffices to show 7 Zc G 53(7? d ) <g> 53(/C) for any closed 
ball C. For this, remember that IC D := {7L € /C: K fl D = 0} € 53 (/C) for 
any open ball Z7 [see Matheron (1975)]. Thus, the proof of Proposition 1 is 
completed by noting that 

7Z C C = {(x,K) £ R d x /C:K n (C - x) = 0} = (J (B x JC D{B ’ C) ), 

BeB 

where D(B, C) = {u — v :u G C, u G 77} is an open ball for B G £>. □ 

Let there be given an unmarked point process 'F = X)?:>i on R d and 
a random compact set So with distribution Q = P o Eg 1 on [/C,53(/C)]. The 
corresponding independently marked point process 'Fq = X)i>i <5 [a';,h,] on R d 
with mark distribution Q is then defined [see Daley and Vere-Jones (1988) or 
Stoyan, Kendall and Mecke (1995)] to be a random element (over [Q,2l, P]) 
that takes values in [M/c , 51fo] > the distribution of which is uniquely deter¬ 
mined by its probability generating functional Gq, Q H = E(n i >i^,3 i )) = 
G^[vq\ (see Section 2), where vq(-) := J lc v(-,K)Q(dK ) and the function 
v: R d x /C i—» [0,1] is Borel measurable such that 1 -v(-,K) has bounded sup¬ 
port for all K G 1C. In the special case of a stationary independently marked 
Poisson process II^q = Z)j>i the shape of Gn A [-] implies that 

(A.2) 6tt A:( 3 M =exp|A J d J (v(x,K) - l)Q(dK)dx^ 

and, furthermore, the Campbell-type formula 

( A/ \ /c 

E* [ f j (x,K)Q(dK)dx 

J j=l JR JK 

holds for any measurable functions /i,..., /*,: R d x /C i—> [0, oo], where the 
sum Y* on the left-hand side of (A.3) stretches over fc-tuples of pairwise 
distinct indices. 

As announced, we conclude the Appendix by showing that, under the 
assumption E|Eo| < oo, the condition E|So + B e {o)\ < oo for some e > 0 is 
not only sufficient as shown by Heinrich (1992), but even necessary for the 
closedness of the stationary PGM S = E(Ha,q). 

Proposition 2. Let Eo be a compact typical grain of the PGM (1.1) that 
satisfies E|So| < oo and E|Eo + 7? £ (o)| = oo for any e > 0. Then P(S is closed 
in R d ) = 0. 
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Proof. Choose K n £ /C, n > 1, such that K n | as n —+ oo and let 
E n ('0) be defined as in the proof of Proposition 1. Obviously, {ijj £ M/c : E n (?/>)n 
2? e (o) = 0 } |{'i/’£ M/c :E(?/>) O -B e (o) = 0 } as n —> oo. Furthermore, since 
(Ej + Xj) D R e (o) 7 ^ 0 iff — Xi £ Sj + B £ (o), we find, using (A.2), that 

p(H n (n AiQ ) n B e {o) = 0 ) = Ena - Iav^WW-^)) 

i>l 

= exp{— AE|(— AT n ) n (E 0 + B £ (o))\}. 

By the monotone convergence theorem and our assumptions, 
lim E\(-K n ) O (S 0 + -B e (o))| = E|E 0 + .B e (o)| =oo, 

n—>oo 

which means for the stationary PGM E = S(II A n) that 

P(E O B e (o) / 0 ) = 1 - lim P(S n (n A q) O B £ (o) = 0 ) = 1 for any e > 0. 

n — kx> 

Thus, 

P(E is closed) = P | E is closed, Q {E 0 Bi/ m (o) / 0 } 

V m> 1 

< P(o £ E) = 1 — exp{—AE|E 0 |} < 1. 

Similarly, P(S is closed) < P(xi £ S,... ,x n £ 5) for any x \,..., x n £ R d . In 
view of (1.4)—(1.6), the probability p£\xi,...,x n ) is arbitrarily close to 
(P(o £ E)) n whenever the distances between the points xi,... ,x n are suffi¬ 
ciently large. This proves the assertion of Proposition 2. □ 

Acknowledgment. I thank one of the referees for helpful comments and 
suggestions. 
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